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Using the local potential approximation of the exact renormalization group (RG) equa- 
tion, we show the various domains of values of the parameters of the 0(l)-symmetric scalar 
Hamiltonian. In three dimensions, in addition to the usual critical surface S c (attraction do- 
main of the Wilson-Fisher fixed point), we explicitly show the existence of a first-order phase 
transition domain Sf separated from S c by the tricritical surface St (attraction domain of 
the Gaussian fixed point). Sf and S c are two distinct domains of repulsion for the Gaussian 

■ fixed point, but Sf is not the basin of attraction of a fixed point. Sf is characterized by an 
endless renormalized trajectory lying entirely in the domain of negative values of the ip 4 - 

■^j- ■ coupling. This renormalized trajectory exists also in four dimensions making the Gaussian 

| fixed point ultra-violet stable (and the <p\ renormalized field theory asymptotically free but 

with a wrong sign of the perfect action). We also show that very retarded classical-to-Ising 
crossover may exist in three dimensions (in fact below four dimensions). This could be an 

■ explanation of the unexpected classical critical behavior observed in some ionic systems. 

PACS 05.10.Cc, 05.70.Jk, ll.10.Hi 

1 Introduction 

• The object of this paper is to carry on studying the local potential approximation of the exact 
renormalization group (RG) equation for the scalar theory Q. In a previous publication [Q] (to be 
considered as the part I of the present work) , we had already considered this approximation with 
a view to qualitatively discuss the connection between the standard perturbative renormalization 
of field theory (as it can be found in most textbooks on field theory, see for example ||) and 
the modern view jl in which the renormalized parameters of a field theory are introduced as 
the "relevant" directions of a fixed point (FP) of a RG transform. Actually the local potential 
approximation, which allows us to consider all the powers of the field ip on the same footing, is 
an excellent textbook example of the way infinitely many degrees of freedom are accounted for 
in (nonperturbative) RG theory. Almost all the characteristics of the RG theory are involved 
in this approximation. The only lacking features are related to phenomena highly correlated to 
the non local parts neglected in the approximation and when the critical exponent rj is small 
(especially for d = 4 and d = 3), one expects the approximation to be qualitatively correct on 
all aspects of the RG theory ^ . 

In the following we look at the domains of attraction or of repulsion of fixed points in the 
0(1) scalar theory in three and four dimensions {d = 3 and d = 4). 

At first sight, one could think that the issue considered is very simple since, with regard to 
criticality, the 0(l)-symmetric systems in three dimensions are known to belong to the same 
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class of universality (the Ising class). Now, because the Ising class is associated to the domain 
of attraction of the unique (non-trivial) Wilson-Fisher fixed point ||, then by adjusting one 
parameter (in order to reach the critical temperature^]) any O(l) scalar Hamiltonian should be 
driven to the Wilson-Fisher fixed point under the action of renormalization. Consequently there 
would have only two domains for the O(l) scalar theory: the critical subspace S c (of codimension 
1) in the Wilson space (S) of infinite dimensions of the Hamiltonian parameters (in which the 
RG transforms generate flows) and the complement to S of S c (corresponding to noncritical 
Hamiltonians). 

In fact, that is not correct because there is another fixed point in S: the Gaussian fixed point 
which, although trivial, controls tricritical behaviors in three dimensions. Now each FP has its 
own basin of attraction in S |j. The attraction domain of the Gaussian FP is the tricritical 
subspace St of codimension 2 (with no intersection with S c ). In addition, we show that there is 
a second subspace of codimension 1 in S, called Sf, which is different from S c , and thus which 
is not a domain of attraction to the Wilson-Fisher fixed point. There is no FP to which a point 
of Sf is attracted to. Sf is characterized by a negative sign of the (/3 4 -Hamiltonian parameter U4 
and is associated with systems undergoing a first-order phase transition ||. We show that an 
endless attractive RG trajectory is associated to this domain of first-order transitions. It is a 
renormalized trajectory (denoted below by Tj) that emanates from the Gaussian fixed point. 
The frontier between Sf and S c corresponds to the tricritical subspace St which is the domain 
of attraction of the Gaussian fixed point while Sf and S c are two distinct domains of repulsion 
for the Gaussian fixed point. 

Actually, the situation is conform to the usual view. Considering the famous ip -model 
[Landau-Ginzburg- Wilson (LGW) Hamiltonian] in which the associated coupling U4 is positive, 
then the Hamiltonian at criticality is attracted exclusively to the Wilson-Fisher fixed point, but 
if U4 is negative, a (/7 6 -term is required for stability, but then one may get either a tricritical 
phase transition or a second- or a first-order transition J?]]. In the present study we do not 
truncate the Hamiltonian which involves all the powers of the field (p. 

We explicitly show that a system which would correspond to an initial point lying very close 
to the frontier St in the critical side (in S c ) would display a retarded classical-to- Ising crossover 
|| . This result is interesting with regard to ionic systems (for example) in which a classical 
behavior has been observed while an Ising like critical behavior was expected. The eventuality 
of a retarded crossover from the classical to the Ising behavior has previously been mentioned 
but without theoretical explanation on how this kind of crossover could develop ||. In H a 
calculation suggests that the RPM model for ionic systems would specifically correspond to 
a scalar Hamiltonian with a negative sign for the </3 4 -Hamiltonian parameter (but the order 
parameter chosen is not the bulk density [JToJ] ) . That calculation has motivated the present 
study. 

We also indicate that the renormalized trajectory T 1 still exists in four dimensions. This 
makes the Gaussian fixed point ultraviolet stable and the scalar field theory formally asymptot- 



ically free. However the associated "perfect" action [TT|] would have the wrong sign to provide 
us with an acceptable (well defined) field theory. 

The organisation of the paper is as follows. In section Q we briefly present the local potential 
approximation of the exact RG equation to be studied. We introduce the strategy we have 
chosen to solve the resulting nonlinear differential equation with a view to show the trajectories of 
interest in the space S of infinite dimension. Because the practical approach to the Gaussian fixed 
point is made difficult due to the logarithmic slowness characteristic of a marginally irrelevant 
direction (for d = 3), we found it useful to first test our numerical method with a close approach 
to the Wilson-Fisher fixed point. We present the characteristic results of this approach and the 



We assume that the second relevant field, corresponding to the magnetic field for magnetic systems, is set 
equal to zero. 
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various kinds of domains corresponding to 114 > (somewhat a summary of Q) 

In section |3] we describe the various kinds of attraction or repulsion domains of the Gaussian 
fixed point (for a negative value of the </? 4 -Hamiltonian parameter) corresponding to tricritical, 
critical and first-order subspaces. Then we discuss the consequences and especially explicitly 
show how a retarded crossover from the classical to the Ising behavior can be obtained. 
We then shortly discuss the case d = 4 when U4 < 0. 

In two appendices we report on some technical aspects of the numerical treatment of the RG 
equation studied, in particular on the appearance of spurious nontrivial tricritical fixed points 
(Appendix |A.1|) . 



2 The RG equation studied 



The local potential approximation has been first considered by Nicoll et al [12| from the sharp 



cutoff version of the exact RG equation of Wegner and Houghton [13], it has been rederived 
by Tokar [14] by using approximate functional integrations and rediscovered by Hasenfratz and 



Hasenfratz [15]. As in (2| we adopt the notation of the latter authors and consider the following 
nonlinear differential equation for the simple function f(l,(p): 



f 



K d f" 
2 1 + /' 



+ 1 



(i) 



in which a prime refers to a derivative with respect to the constant dimensionless field 99 (at 
constant I) and f(l,(f) = V'(l,tp) is the derivative of the dimensionless potential V(l, (p); f 
stands for df /dl\ v in which I is the scale parameter defined by =e~ l and corresponding to 
the reduction to A of an arbitrary initial momentum scale of reference Ao (the initial sharp 
momentum cutoff). Finally, is the surface of the d-dimensional unit sphere divided by (27r) rf . 

A fixed point is solution of the equation / = 0. The study of the resulting second order 
differential equation provides the following results: 



d > 4, no FP is found except the Gaussian fixed point. 

3 < d < 4, one nontrivial FP (the Wilson-Fisher fixed point ||) is found [^, 16, |l7]] 



A new nontrivial FP emanates from the origin (the Gaussian fixed point) below each 
dimensional threshold dk = 2k/(k — 1), k = 2, 3, . . . , 00 fL8f. 



If one represents the function f(l,<p) as a sum of monomials of the form: 



n 



then, for d = 3, the Wilson-Fisher fixed point /* is located in S at ||: u% = — 0.461533 •• •, 
u\ = 3.27039 ■ • • , u% = 14.4005 ■••,«! = 32.31289 • ■ ■ , etc. 

Once the FP is known, one may study its vicinity which is characterized by orthogonal 
directions corresponding to the infinite set of eigenvectors, solutions of the differential equation 
(jl|) linearized at /*. The eigenvectors associated to positive eigenvalues are said relevant; when 



the eigenvalues are negative they are said irrelevant and marginal otherwise [19|. 

The relevant eigenvectors correspond to directions along which the RG trajectories go away 
from the FP and the irrelevant eigenvectors correspond to directions along which the trajectories 
go into the FP. A marginal eigenvector may be relevant or irrelevant. 

Our present FP /* has only one relevant direction and infinitely many irrelevant directions 
(no marginal direction, however see BJ). As already explained and shown in pf, in order to 
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approach /* starting from an initial point in S, one must adjust one parameter of the initial 
function /(0, (f). This amounts to fixing the temperature of a system to its critical temperature. 
Starting with a known initial function (at "time" I = 0) say: 

f(0,<p) =U 2 (0)^ + U 4 (0)V? 3 

we adjust U2 (0) to the critical value = —0.29958691 • • • corresponding to (0) = 3 so that 
f(l,tp) [solution at time I of the differential equation ([]])] approaches /* when / — > oo. The 
approach to /* is characterized by the least negative eigenvalue A2 = —1/loi {oj\ was noted to in 
IH). This means that, in the vicinity of /* any parameter u n (I) evolves as follows (/ — ► 00): 

u n (I) ~ u* n + a n exp {—u)il) 

Fig. H illustrates this feature for the first four u n (Z)'s in the approach to /*. In Q the 
two associated attractive trajectories (locally tangent to the least irrelevant eigenvector in the 
vicinity of /*) was noted Ti and T[. 

One may also constrain the trajectory to approach /* along the second irrelevant direction 
(with the associated attractive trajectories noted T2 or T 2 in Q and associated with the second 
least negative eigenvalue A3 = — l/o^)- I n this case a second parameter of the initial / must be 
adjusted, e.g., U4 (0) must be adjusted to u\ and simultaneously 112 (0) to the corresponding uf, 
see [jl6] , [2|. Then, in the vicinity of /*, any parameter u n (I) will evolve as follows: 

u n (I) =u* n + a' n exp {-UJ2I) 

Looking for this kind of approach to /*, we have found that 6.66151663 < u\ < 6.66151669 
and u% = -0.58328898880579 • • • This has allowed us to estimate uj 2 = 2.84. Although the 
shooting method is certainly not well adapted to the determination of the eigenvalues (see the 
huge number of digits required in the determination of u\ and u%), our estimate is close to 



UJ2 — 2.8384 found by Cornelias and Travesset [20]. 



Because u\ cannot be perfectly determined, the trajectory leaves the trajectory T2 before 
reaching /* to take one of the two directions Ti or T 1 (corresponding to Fig || illustrates 
this effect with the evolution, for n = 2, of the following effective eigenvalue: 

in) m _ d 2 u n (l)/dl 2 
^ {l) - " du n (l)/dl (2) 
the definition of which does not refer explicitly to /*. The evolution of uj c q(1) shows a flat 
extremum (or a flat inflexion point) at an RG eigenvalue of /* each time the RG flow runs along 
an eigendirection in the vicinity of /*. 

Similarly to u%, the value u\ cannot be perfectly determined, consequently the trajectory ends 
up going away from the fixed point. This provides us with the opportunity of determining the 
only positive (the relevant) eigenvalue corresponding to the critical exponent v = \\ = —1/ujq 
[uj e g(l) shows then a flat extremum at loq when the flow still runs in the close vicinity of /*]. 
Finally, far away from the fixed point, the RG trajectory approaches the trivial high temperature 
fixed point characterized by a classical eigenvalue (equal to i). The global picture summarizing 
the evolution of cj e ff (0 along a RG trajectory initialized in such a way as to approach /* first 
along T2, is drawn on fig. ||[ 

The values we have determined by this shooting method are (for eigenvalues other than u>2 
already mentioned): 

wi = 0.5953 
v = 0.68966 



which are close to the values found, for example, in |lj|, gfj: u\ = 0.5952 and v = 0.6895 
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3 Trajectories for w 4 < 



In the preceding section, we have obtained a RG trajectory approaching the Wilson-Fisher fixed 
point /* along T2 by adjusting two parameters of the initial Hamiltonian (u% and Uq). This 
is exactly the procedure one must follow to determine a tricritical RG trajectory approaching 
the Gaussian fixed point in three dimensions (because of its two relevant directions). The only 
difficulty is to discover initial points in S which are attracted to the Gaussian fixed point. To 
this end, we again use the shooting method. 

From usual arguments on the LGW Hamiltonian and from the work done by Aharony on 
compressible ferromagnets ||, one expects to find the tricritical surface in the sector 114 < 
(and with 112 > 0). Thus we have tried to approach the Gaussian fixed point starting with initial 
function /(0) of the form: 

/(0, cp) = u 2 (0) <p + u 4 (0) <p 3 + u 6 (0)</ (3) 

with (not large) negative values of u 4 (0), for example U4 (0) = —1. 

Because the Gaussian fixed point is twice unstable, we must adjust two parameters to ap- 
proach it starting with (H|). We do that by successive tries (shooting method). For example, if 
we choose 114 (0) = — 1 and uq(0) = 3 and determine a value of 1*2 (0) such as to get a trajectory 
which does not go immediately towards the trivial high temperature fixed point, the best we 
obtain is a trajectory which approaches the Wilson-Fisher fixed point (thus the corresponding 
initial point belongs to the attraction domain of /* although 114 (0) < H). But if uq(0) = 2, 
the adjustment of 112 (0) with a view to counterbalance the effect of the most relevant direction 
of the Gaussian fixed point (which would drive the trajectory toward the high temperature FP) 
yields a runaway RG flow towards larger and larger negative values of U4 (I). From now on, 
the target is bracketted: the tricritical trajectory corresponding to U4 (0) = — 1 can be obtained 
with a value of Ug(0) in the range ]2, 3[ (we actually find a rather close approach to the Gaussian 
fixed point for 2.462280 > i* 6 (0) > 2.4622788 and 6.4618440- • • > 1*2(0) > 6.4618407- • •). 

In order to understand the origin of the direction of runaway in the sector of negative values 
of U4, it is worth to study the properties of the Gaussian fixed point by linearization of the RG 
flow equation in the vicinity of the origin. If we request the effective potential to be bounded by 
polynomials then the linearization of eq. ((!]) identifies with the differential equation of Hermite's 



polynomials of degree n = 2k — 1 for the set of discrete values of A satisfying [15 

2 + d - 2A 



d-2 

from which it follows that 



2k -I fc = l,2,3,--- (4) 



• for d = 4: Afc = 4 — 2k k = 1, 2, 3, • • there are two non-negative eigenvalues: Ai = 2 
and A2 = 

• for d = 3: A& = 3 — k k = 1, 2, 3, ■ • -, there are three non-negative eigenvalues: Ai = 2, 
A 2 = 1 et A 3 = 

If we denote by Xki}f) the eigenfunctions associated to the eigenvalue Afc, it comes: 

• Xi = xt = V? 3 — I'/ 9 ' X3 = ~ + x^ 3 ' ' ' '' whatever the spatial dimensionality d. 

The upperscript "+" is just a reminder of the fact that the eigenfunctions are defined up to 
a global factor and thus the functions Xki.^) = ~Xk( l P) are a ^ so eigenfunctions with the same 
eigenvalue A&. 
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3.1 Case d = 3 



Similarly to xt i the direction provided by X2 in S 1 is a direction of instability of the Gaussian 
fixed point. Now xt ls associated with the well known renormalized trajectory Ti on which is 
defined the usual (massless) 993-field theory Ji~6|, |2|, for the same reasons a renormalized trajectory 
T 1 locally tangent to X2 i n the vicinity of the origin of S exists with the same properties as Tj 
(see Q). The difference is that T 1 lies entirely in the sector U4 < and is endless (not ended 
by a fixed point). 

This endless renormalized trajectory is associated with systems undergoing a first-order phase 
transition. This is due to the absence of fixed point [|21jJ, in which case the correlation length 
£ cannot be made infinite although for some systems lying close to T x and attracted to it (i.e. 
at the transition temperature), £ may be very large (because T 1 is endless), in which cases one 
may say that the transition is almost of second order p2| . Of course, a domain of first-order 
phase transition in S was expected from the usual arguments ]?], 0], we only specify better the 
conditions of realization, in S, of the first-order transition. 

Fig. U shows the attractive trajectory T 1 together with the attractive tricritical line ap- 
proaching the Gaussian fixed point. The tricritical surface St separates the first-order surface Sf 
from the critical surface S c . Fig. |] shows also that systems lying close to the tricritical surface 
may still be attracted to the Wilson-Fisher fixed point. In this case the effective exponents may 
undergo a very retarded crossover to the asymptotic Ising values compared to usual systems 
corresponding to initial points chosen in the sector 114 > of S. Fig. [5] illustrates how minus the 
inverse of (||) provides us with different evolutions [calculated from ([[])] of the effective exponent 
v es( T ) [with t ex (T — T c )/T c ] according to the initial point chosen in S. It is worth to explain 
how we have defined v e s(r). 

We have seen at the end of section [2] that the quantity (^) undergoes a flat extremum (or a 
flat inflexion point) at an RG eigenvalue of /* each time the RG flow runs along an eigendirection 
in the vicinity of /*. Now it happens that this extremum is less and less flat as one chooses 
larger and larger values of (^2(0) — u^) (for the eigenvalue u) but still exists. This provides us 
with a way to express the evolution of an effective exponent v e $ when the RG-substitute to t, 
namely (112(0) — u^/u^, is varied. Fig. ^ shows such an evolution for some initial Hamiltonian 
(with 114(0) = 4). Notice that for such a Hamiltonian, the extremum disappears before v e s 
reaches the trivial value | (associated with the approach to the trivial high temperature fixed 
point and to a regular — non critical — behavior) while in the case of a Hamiltonian initialized 
close to the tricritical surface, the classical-to-Ising crossover is complete (see fig. ^|). This is 
because in the latter case the RG trajectory comes close to the Gaussian fixed point (and v e &(r) 
has an extremum at ^) before approaching to /*. This reinforces the idea that the so-called 
classical-to-Ising crossover actually exists only between the Gaussian and Wilson-Fisher fixed 
points [f23| . 

The same configuration displayed by fig. || has been obtained also by Tetradis and Litim [ 24 1 



while studying analytical solutions of an exact RG equation in the local potential approximation 
for the 0(iV)-symmetric scalar theory in the large N limit. But they were not able to determine 
"the region in parameter space which results in hrst-order transitions" [H| . 

3.2 Case d = 4 

To decide whether the marginal operator (associated with the eigenvalue equal to zero, i.e. A2 
in four dimensions, or A3 in three dimensions) is relevant or irrelevant, one must go beyond the 



linear approximation. The analysis is presented in [15| for d = 4. If one considers a RG flow 
along xt such that g2(f,l) = c(l)x%( i P)i then one obtains, for small c: c(l) = c(0) [1 — Ac(0)l} 
with A > 0. Hence the marginal parameter decreases as I grows. As it is well known, in four 
dimensions the marginal parameter is irrelevant. However, if one considers the direction opposite 
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to xt (i- e - X2 ) then the evolution corresponds to changing c — ► — c. This gives, for small values 
of c: c(l) = c(0) [1 + A |c(0)| Z] and the parameter becomes relevant. The parameter c is the 
renormalized 4 coupling constant and it is known that in four dimensions the Gaussian 
fixed point is IR stable for ur > but IR unstable for ur < |25[] . 

We have verified that the trajectory T x survives when d = 4 (contrary to Ti, see Q). That 
trajectory T 1 is a renormalized trajectory on which we could define a continuum limit for the 
V?|-field theory and if the corresponding (perfect) action was positive for all 92, one could say that 
the 1^4-field theory with a negative coupling is asymptotically free. Unfortunately, because the 
y? 4 -term is dominant for large tp in the vicinity of the origin of S (due to the relevant direction 
provided by X2O, the negative sign of the renormalized coupling prevents the path integral to 
be well defined. 

However, because the action to which one refers in the continuum limit (the perfect action) 
is formal (because it involves an infinite number of parameters and cannot be written down, 
see §|) we wonder whether the wrong sign of the action is actually a valid argument to reject 
the c/?4-field theory with a negative renormalized coupling. It is worth to mention that the 
asymptotically free scalar field theory which has recently been considered on a lattice [26] could 
actually be the <?!> 4 -field theory with a negative coupling to which we refer here. 

Acknowledgments We dedicate this article to Professor Yukhnovskii in grateful recognition 
of his efficient and generous help in fostering the Ukrainian-French Symposium held in Lviv in 
february 1993, with the hope that in the future the contacts between our two communities will 
further develop. 
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A The finite difference method used 



For technical reasons, instead of studying Eq. (||), we consider the differential equation satisfied 
by g(<p) = f'(f) (i.e. the second derivative of the potential with respect to the field): 



Kd 
2 



g" (g') 



.1 + 5 (i + g)\ 

Starting with a known initial function (at "time" I = 0), we follow its evolution in S by 
approximating the differential equation (|5|) by finite differences and a two dimensional grid with 
the uniform spacings dy = 0.01 and dl = 0.000390625. The finite difference formulas for the 
derivatives g" and g' have been chosen with the accuracy 0(dy 4 ): 



g'(y) = -^[g(y + dy)-g(y-dy)-g(y + 2dy)+g(y-2dy)] + 0(dy 4 ) (6) 
16 

9 (y) = j^2[g{y + dy)+g(y-dy)-30g(y) 

-g(y + 2dy) - g(y - 2dy)} + 0(dy 4 ) (7) 

The evolutionary function g(y,l) is known (calculated) at the discret set of points yi = i ■ dy 
with (i > 0) and a maximum value i max = 82. (This value is large enough to study the approach 
to the Wilson-Fisher fixed point with a great accuracy but is too small for studying precisely 
the approach to the Gaussian fixed point.) At each time If. = k-dl, the derivatives are estimated 
from g(y u l k ) = g(y h lk-i)+g(yi, h-l)-dl by using Eqs(j|,0) which apply only for 1 < i < i max -i. 
For the marginal points i = 0, 1 we use the parity of g(y) [by inserting g(—n ■ dy) = g(n ■ dy) 
for n = 1, 2 in Eqs(^J^)]. For the two other marginal points i = i m iw-i, W of the grid, there is 
no fixed solution and we shall used alternately the two following conditions [using the obvious 
abreviation g(i) instead of g(i ■ dy)]: 



g\i) = g'(i-l) 
g"(i) = g"(i-l) 



g'(i) = i 



1 

dy 



25 4 
—g(i) - 4g(i - 1) + 3g(i - 2) - -g(i - 3) 



1 



4) + 0(dy 4 



9"H) = i 



1 

dy 2 

61 , 5 , 



915 ,., 77 ,. 107 ,. 

q(l) qli — 1) H q(i 

244 yV ; 6 yV J 6 yy 



+ 0(dy 4 



2) - 13g{i - 3) 



(8) 



0) 



Condition 2 is more accurate than condition 1 but leads sometimes to strong unstabilities 
which do not appear when we first use Condition 1 and then Condition 2 after some finite "time" 
Iq. The validity of the procedure is tested by, for example, trying to approach a given fixed point 
(see the main part of the paper). 
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A.l Appearance of spurious fixed points in the approach to the Gausian fixed 
point 

In trying to determine the attractive tricritical trajectory (approaching the Gaussian fixed point), 
we have encountered a spurious twice unstable fixed point lying at some finite and non negligible 
distance to the Gaussian fixed point. To understand the origin of this undesirable numerical 
effect, it is necessary to discuss a bit the solution of the fixed point equation / = 0. 

From flU) or (|5|), one sees that the fixed point equation is a second order non linear differential 
equation and a solution would be parametrized by two arbitrary constants. One of these two 
constants may easily be determined: since g*{<^) is expected to be an even function of <p [0(1) 
symmetry] then g*'(0) = may be imposed. It remains one free parameter, thus a one-parameter 
family of (nontrivial) fixed points are solutions to the differential equation. But there is not an 
infinity of physically acceptable fixed points; all but a finite number of the solutions in the family 
are singular at some ip c |?]]. Formally, by requiring the physical fixed point to be defined 

for all (p then the acceptable fixed points are limited to the Gaussian fixed point and (for d = 3) 
to the Wilson-Fisher fixed point. 

However, in our study, because we numerically consider the function g((f) in some finite 
range of values of ip (see above: i max = 82), it appears that in approaching the origin of S, 
infinitely many pseudo-fixed points exist which have there (^-singularity located outside the 
finite range explicitly considered and there is at least one of them which looks like a tricritical 
fixed point. When we enlarge the range of <p to i max = 200, the previously observed nontrivial 
tricritical fixed point disappeared to the benefit of another one located closer to the origin. In 
conclusion, a larger and larger number of grid-points must be considered as one tries to come 
closer and closer to the Gaussian fixed point. This particularity together with the slowness of 
the approach along a marginal direction makes it excessively difficult to come very close to the 
Gaussian fixed point. 
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Figure captions 



1. Evolutions for d = 3 of the first four Hamiltonian parameters U2(l), 114(1), uq(1), ug(l) in a 
close approach to the Wilson-Fisher fixed point /* along Tj or T 1 . The effective inverse 
eigenvalue u; e ff (I) is given by eq. @ for n = 2, 4, 6, 8. All these quantities reach the same 
universal value uj\ characteristic of the least irrelevant eigendirection of /*. To get this 
close approach to /* from eq. (||), the initial critical value u% corresponding to 114(0) = 3, 
has been determined with more than twenty digits. 

2. When a second condition is imposed on the initial Hamiltonian parameters, the approach 
to /* may be adjusted such as to asymptotically take the second least irrelevant eigendi- 
rection. Here u e g (I) is given by eq. (|2|) for n = 2 it clearly undergoes (full line) a flat 
inflexion point at the value LO2 = 2.84 corresponding to an approach to /* along T2, the 
greater the critical parameter u\ is accurately determined, the longer is the flat extremum. 
Because u\ is not completely determined [within the available accuracy in solving eq. ([l])] 
the trajectory leaves the direction of T2 to take one of the two directions of approach 
associated to the least irrelevant inverse eigenvalue lo\ (corresponding to Ti or T 1 as in- 
dicated by dashed curves). Here, the trajectory corresponding to the full line goes along 
T 1 . Again a flat externum of ui e s (I) indicates the approach along an eigenvector of /* and 
requires an accurate determination of the critical value u§. Because this determination 
is not complete, the trajectory ends up going away from /* as indicated by the sudden 
departure of u; e ff (I) from uj\ for the large values of I. 

3. This figure is a continuation of fig. ||. It shows the various plateaux that u) e g (I) undergoes 
along a RG trajectory first adjusted to approach /* along the second irrelevant direction 
(plateau at UJ2 = 2.84). Because it is not possible to determine exactly the initial condi- 
tions, the trajectory always ends up going away the fixed point towards the trivial high 
temperature fixed point characterized by the classical value ^ (for minus the inverse of 
u! e s (I), thus the final plateau at —2). In-between, the RG flow has been influenced by the 
close vicinity of the least irrelevant eigenvector (plateau at U\) and that of the relevant 
eigenvector (plateau at ujq = — jj). The various regimes of the RG flows are indicated by 
the vertical arrows on the left (direction of the flow with respect to the fixed point) and 
on the right of the figure (distance to the fixed point). 

4. Domains of attraction and repulsion of the Gaussian fixed point. The figure represents 
projections onto the plane {^21^4} of various RG trajectories running in the space S 
minus one dimension. The flows have been obtained by solving Eq. (|l|). Black circles 
represent the Gaussian and the Wilson-Fisher (W-F FP) fixed points. The arrows indicate 
the directions of the RG flows on the trajectories. The ideal trajectory (dot line) which 
interpolates between these two fixed points represents the usual renormalized trajectory 
Ti corresponding to the so-called 993 renormalized field theory in three dimensions (usual 
RT for U4 > 0). White circles represent the projections onto the plane of initial critical 
Hamiltonians. For 7/4(0) > 0, the effective Hamiltonians run toward the Ising fixed point 
asymptotically along Tj (simple fluid). Instead, for 1*4(0) < and according to the values 
of Hamiltonian coefficients of higher order (uq, uq, etc.), the RG trajectories either (A) 
meet an endless RT emerging from the Gaussian FP (dashed curve) and lying entirely 
in the sector M4 < or (B) meet the usual RT Ti to reach the Ising fixed point. The 
frontier which separates these two very different cases (A and B) corresponds to initial 
Hamiltonians lying on the tri-critical subspace St (white square C). This is a source of RG 
trajectories flowing asymptotically toward the Gaussian FPalong the tricritical RT. Notice 
that the coincidence of the initial point B with the RG trajectory starting at point A is 
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not real (it is accidental, due to a projection onto a plane of trajectories lying in a space of 
infinite dimension). The points A or B could correspond to the restricted primitive model 
of ionic systems (see Q). 

5. Evolutions of an effective exponent v e fi( r ) [with r oc (T — T c )/T c ] along three different 
families of RG trajectories (see text for additional details). The full squares indicate the 
evolution of v e s(r) for a family of trajectories initialized in the sector U4 > with 1*4(0) = 3 
and for various values of 1*2(0) (the same system at criticality corresponds to the white 
circle "Simple fluid" of fig. Q). When 1*2(0) — » u\ the effective exponent approaches the 
critical exponent value v = 0.69 compatible with the present study. One observes that 
the crossover towards the classical value ^ is not complete because v e fi(r) ceases to make 
sense before r becomes large. This is not the case of the evolution represented by the full 
circles which corresponds to trajectories initialized close to the Gaussian fixed point. In 
this case the complete crossover reproduces the interpolation between the Gaussian and 
the Wilson-Fisher fixed points and typically corresponds to the usual answer given by field 
theory |28[j . The third evolution (full triangles) corresponds to a family of Hamiltonian 
initialized close to the tricritical surface but still attracted to the Wilson-Fisher fixed point. 
One sees that the classical-to-Ising crossover is complete but highly retarded compared to 
the two other cases. This is because at criticality, the RG flow is first attracted to the 
Gaussian fixed point (showing then an apparent classical value of v) before interpolating 
between the Gaussian and the Wilson-Fisher fixed point. 

6. Illustration of the evolution of the extrema u e g(l) [minus the inverse of eq. @] for various 
values of r = (1*2(0) — u|) ju\ and for the family of RG flows initialized at 1*4(0) = 3. The 
extremum (grey triangle) disappears at some not very large value of r (about 10 -0 ' 5 ) and 
does not reach the classical value 4. This induces the partial Ising-to-classical crossover 
drawn on fig. || (squares). 
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